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We analyze spin scattering in ballistic transport of electrons through a ripple at a normal incidence 
of an electron flow. The model of a ripple consists of a curved graphene surface in the form of 
an arc of a circle connected from the left-hand and right-hand sides to two flat graphene sheets. 

At certain conditions the curvature induced spin-orbit coupling creates a transparent window for 
incoming electrons with one spin polarization simultaneously with a backscattering of those with 
opposite polarization. This window is equally likely transparent for electrons with spin up and spin 
down that move in opposite directions. The spin filtering effect being small in one ripple becomes 
prominent with the increase of N consequently connected ripples that create a graphene sheet of the 
sinusoidal type. We present the analytical expressions for spin up (down) transmission probabilities 
as a function of N connected ripples. 

PACS numbers: 72.25.-b,71.70.Ej,73.23.Ad 


I. INTRODUCTION 

The extraordinary properties of graphene have at¬ 
tracted enormous experimental and theoretical attention 
for a decade (see e.g. RefsIUi). Graphene being a zero- 
gap semiconductor has a band structure described by a 
linear dispersion relation at low energy, similar to mass¬ 
less Dirac-Weyl fermions. Such a band structure leads to 
exceptionally high mobility of charged carriers. A ques¬ 
tion of possible mechanisms that would allow us to throt¬ 
tle the mobility and, consequently, to control a conduc¬ 
tivity is a topical subject in graphene physics, due to its 
fundamental as well as technological significance. 

Among various mechanisms that might affect the mo¬ 
bility, the scattering that could be induced by a ripple 
(see, for example, discussion in RefH) appears to be the 
most natural one, since graphene sheets are not perfectly 
flat. Moreover, periodic ripples can be created and con¬ 
trolled in suspended graphene, in particular, by thermal 
treatment^ and by placing graphene in a especially pre¬ 
pared substrate. Indeed, curvature of the surface affects 
the TT orbitals that determine the electronic properties of 
graphene. It results in enhancement of spin-orbit cou¬ 
pling that could serve as a source of spin scattering. We 
recall that the intrinsic (intraatomic) spin-orbit interac¬ 
tion in flat graphene is weabi*^!^. It makes spin deco¬ 
herence in such a material weak as well, i.e., scattering 
due to disorder is supposed to be unimportant. In order 
to get deep insight into the nature of curvature induced 
scattering, it is desirable to elucidate among many ques¬ 
tions the basic one: What are the distinctive features of 
curvature induced spin-orbit coupling ? One can further 
ask how to employ these features to guide an electron 
transport in a graphene-based system at the theoretical, 
and, quite likely, practical levels. 

A consistent approach to introduce the curvature 
induced spin-orbit coupling (SOG) in the low energy 


physics of graphene have been developed by Ando^ and 
by othersi"— in the framework of effective mass and tight- 
binding approximations. Recent measurements in ul¬ 
tra clean carbon nanotubes (GNTs)^, i.e., in an ex¬ 
treme form of curved graphene, revealed the energy 
splitting that can be associated with spin-orbit cou¬ 
pling. The measured shifts are compatible with theo¬ 
retical predictions^!^, while some features regarding the 
contribution of different spin-orbit terms in metallic and 
non-metallic GNTs are still debatable (see, for example, 
discussion inli7i£). Nowadays, nevertheless, there is a 
consensus that for armchair GNTs one obtains two SOG 
terms: one preserves the spin symmetry (a spin projec¬ 
tion on the CNT symmetry axis), while the second one 
breaks this symmetry^d^d^d^. Thus, we have a reliable 
answer to the first question, at least, for armchair GNTs. 
In some previous studies^i^did^ the role played by the 
second term was underestimated. In this paper we will 
attempt to answer how full curvature induced SOC, in¬ 
cluding the second term, could be used to create a po¬ 
larized spin current with a high efficiency in a rippled 
graphene system. 

The structure of this paper is as follows. In Sec. II 
we briefly discuss the explicit expressions for the eigen 
spectrum and eigenfunctions of an armchair nanotube 
with a full curvature induced spin-orbit coupling. By 
means of these results we introduce a scattering model 
for one ripple and extend this model for N continuously 
connected ripples. In Sec. Ill we provide a discussion 
of our results in terms of simple estimates. The main 
conclusions are summarized in Sec. IV. 


II. SCATTERING PROBLEM 

In order to model a scattering problem on a ripple we 
consider a curved surface in the form of an arc of a circle 
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connected from the left-hand and right-hand sides to two 
flat graphene sheets. The solution for flat graphene is 
well knowi>ii>^. The solution for a curved graphene sur¬ 
face can be expressed in terms of the results obtained for 
armchair CNTs in an effective mass approximation, when 
only the interaction between nearest neighbor atoms is 
taken into accounti^. 


frame, and obtains 


H' = U{e)HU-\e) = Hmu + Hsoc , 

Hkin = -17 \Ty® Idy + Tx® I-^de 


Hsoc = 




V ■ 


( 6 ) 


A. Low energy spectrum of the armchair nanotube 


Here the operators Tx,y,z are the Pauli matrices that act 
on the wave functions of A- and B-sub-lattices (a pseudo¬ 
spin space), and 


Let us recapitulate the major results^^ in the vicinity 
of the Fermi level E = 0 for a point K in the presence of 
the curvature induced spin-orbit interaction in an arm¬ 
chair CNT. The y axis is chosen as the symmetry and the 
quantizations axis. In this case the eigenvalue problem 
is defined as 



with the following definitions: 

f = j{kx-iky) + i^axif)-^-^ay, (2) 

k — i ^ k — —\ — 

(f^) = cos 9 — Gz sin 9 . 


Here, d'x,y,z are standard Pauli matrices, and the spinors 
of two sub-lattices are 


t^K 


rpK 

rpK 

^AA 


F, 


K 


^ BA 
pK 


( 3 ) 


The following notations are used: 7 = —y/3Vppaj2 = 

7 oa, 7 ' =_^{V^p - V^p)al2 = 71 a, p=l- 87787 (see 
e.g. Refla). The quantities and are the trans¬ 
fer integrals for a and tt orbitals, respectively in a flat 
graphene; a = y/Sd ~ 2.46 A is the length of the prim¬ 
itive translation vector, where d is the distance between 
atoms in the unit cell. 

The intrinsic source of the SOC 6 = A/(3e,r(T) is de¬ 
fined as 


. . 3 h dV ^ dV ^ , . 


( 4 ) 


where V is the atomic potential and e-^cr = ejp ~ ^ 2 p- Tii® 
energy ejp is the energy of u-orbitals, localized between 
carbon atoms. The energy ejp is the energy of 7 r-orbitals, 
directed perpendicular to the curved surface. 

By means of the unitary transformation 

U(9) = exp{i^ay) ® I, (5) 

where / is 2 x 2 unity matrix, one removes the 9 depen¬ 
dence in Hamiltonian o, transformed in the intrinsic 


A, = 7 (1 + 4<5p)/(2i?), \y= 6 -t'/{AR) (7) 


are the strengths of the SOC terms. In the Hamiltonian 
(ini) the term A^;) conserves, while the other one Ay) 
breaks the spin symmetry in the armchair CNT. 


The operator Jy, being an integral of motion 


H,Jy 


0 , is defined in the laboratory frame as 


Jy=/®(^Ly + ^)=J0(-i5, + ^) , (8) 

while in the intrinsic frame it is 


Jy^j'y = UJyll-^ = I® {-ids) ■ (9) 

This integral allows to present the wave functions as 


F'{9,y) = 



( 10 ) 


These wavefunctions are also the eigenfunctions of the 
other integral of motion, the operator k'y = ky. Here, 
m = ± 1 / 2 , ±3/2,.... is an eigenvalue of the angular mo¬ 
mentum operator Jy. For the components of the eigen¬ 
vector F'{9,y) the relations |A| = |iA| and \B\ = \C\ are 
fulfilled at real values of m and ky. 

Solving the eigenvalue problem H'F' = EF', one ob¬ 
tains the eigen spectrum 


E = zLEm,q , Ezn^q = + ty + Xy + + 2Dzn,q , 

Dm,q = ^x + A) (H) 

where tm = 'ym/R, ty = ^ky. 


B. Scattering model for one ripple 

Keeping in mind a discussion that will be given here¬ 
after, we analyze the following geometry (see the con¬ 
struction profile on FiglT|). It consists of one arc of a 
circle that is connected from the left-hand side to a flat 
graphene sheet. This (direct) arc is continuously con¬ 
nected to the inverse arc of the same radius that is con¬ 
nected to the right-hand flat graphene sheet. We put the 
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FIG. 1: The rippled graphene system. 


origin of the coordinate at the center of the direct arc of 
the circle. 

To give better insight into the scattering phenomenon 
in our model of a ripple, we study first only the direct 
arc of the circle connected to two flat surfaces. Two 
flat surfaces are: i)the region L, defined in the intervals 
—oo<x< —R cos 00] the region R, defined in the inter¬ 
vals i? cos 00 < X < oo. The region I is a part of a nan¬ 
otube of radius R, defined as — i?cos0o < x < i?cos0o- 
At 00 = 0, the ripple is a half of the nanotube, while at 
^0 = 7r/2 the ripple does not exist. For the sake of anal¬ 
ysis we introduce the angle (p = tt — 20o. To describe the 
scattering phenomenon one has to define wave functions 
in different regions: flat (L,R) and curved (I) graphene 
surfaces. 

Regions L and R are described by the Hamiltonian 

Ho =1 {fxkx + Tyky) ® I (12) 

that does not mix spin components. For the sake of sim¬ 
plicity, we consider the electron motion at the normal in¬ 
cidence, with the electron wave vector k = {kx,0). One 
solves the stationary Schrodinger equation Ho'i’ = Eo'i' 
and obtains the corresponding eigenstates 

Eo = ±^\kx\, (13) 

T' = exp {ikxx)%(k ^), (14) 

N _ 1 sign{jkxEo)% 

“ 2 I 

= ,a = ±l. (16) 

Evidently, the wave functions in regions L, R, can be 
written as a superposition of all possible solutions for flat 
graphene. To proceed, with the aid of eigenspinors (ITSl) . 
m, we introduce an auxiliary matrix Mq (4 x 4) for a 
given value of energy at the normal incidence 



Mo = ,^o\Kx) ,'fo\-K.)) . 

(17) 

Here, we define the variable = sign( 7 Eo) |fca; | to ensure 
that the first two columns of the matrix Mq correspond 
to eigenstates that move in a positive cc-direction, while 
the last two columns correspond to eigenstates that move 
in a negative a:-direction. 



m 


FIG. 2: (Color online) The spectrum (I19II at = 0 as a 
function of the quantum number m. Dashed (red) and solid 
(blue) lines are associated with states characterized by the 
quantum number nis^-i and ms=+i, respectively. The val¬ 
ues of ±ma at the energy Ef = 0.2eV (solid horizontal line) 
are indicated by vertical dashed lines. The parameters of cal¬ 
culations are R — lOA, S = 0 . 01 , p = 0 . 1 , 7 = |l.42AeV, 

7 ' = 7 |, A, = i(l/2 -b 2Sp) = 0.32eV, \y = ^ ^ 0.0043eV. 
The arrow indicates the gap displayed on the insert. 


The matrix Mq is unitary, i.e., = Mq. It allows 

us to define a general form of the wave function for 
a flat graphene 

'i’L.Rix) = Mo exp (^ik{x - xl,r)'^ Cl,r , (18) 

where K = dia.g{Kx, —Kx, —K^) is a diagonal ma¬ 
trix, xl,r are a;-coordinates where flat and curved sur¬ 
faces are connected, and Cl,r are corresponding vectors 
with four unknown yet, normalized coefficients in each 
region. Note that we do not consider inelastic scattering. 
Therefore, since the electron energy is conserved, we use 
the same vector k = {kx,0) for the left and right flat 
graphene surfaces. 

For the curved surface we use eigenspinors of the 
Hamiltonian ([6|). The general form of these eigenspinors 
is defined in the intrinsic frame^i. Therefore, we apply 
the inverse transformation ([5]) to these eigenspinors in or¬ 
der to analyze the scattering problem in the laboratory 
frame. At ky = 0 the spectrum m and the eigenspinors 
are particularly simple 


Ea = + sXx , s = ±l, 

4' = exp(im0) , 

(t ) — f \ 


= exp (-if (Ty) 


itm 

(^Ea Xx) 


(19) 

( 20 ) 
( 21 ) 

( 22 ) 


Note that energies in flat graphene and in a curved sur¬ 
face are different Eo = Ea+A^^)"^ (see details in Reflisll. 
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This effect is caused by different hybridizations of tt elec¬ 
trons in flat graphene and a graphene-based system with 
curvature. In particular, A = 5/6, 7/12 (eV) in the arm¬ 
chair and zig-zag nanotubes, respectively. 

At a fixed energy of the electron flow Eq Ea, 
Eg. lfT^ yields four possible values of the quantum num¬ 
ber m 

R I - 

m ^ TOs = ± —y'(sAa - Aa;)2 - Ag, S = ±1. (23) 

Since the angular momentum is not longer the integral of 
motion, we have to consider the mixture of eigenfunctions 
with all possible values of to at a given energy. 

As an example of the spectrum CH), a few positive 
energy branches as a function of the quantum number 
TO are shown in Fig.2. The branches are distinguished 
by the index s = ±1. There is an anti-crossing effect 
between energy states characterized by the same ms=+i 
quantum number. This anticrossing is brought about by 
the interaction Xy) that breaks the spin symmetry 
(see Sect.IIA) in the curved graphene surface. It results 
in a gap of = 2Xy near Ea = Xx indicated by the ar¬ 
row (see the inset on FigI5]). Similar gap occurs near 
Ea = —Xx for the TO-states with index s = — 1. As a 
consequence of these gaps, evanescent modes arise at en¬ 
ergies Xx — Xy < \E\ < Xx + Xy in our system. For the 
sake of illustration the positive spectrum (IT9t of TO-states 
is crossed by the horizontal line that mimics the Fermi 
energy. The crossing points determine quantum num¬ 
bers TO that have non-quantized values when the curved 
surface (arc of circle) is connected to the flat one. 

With the aid of eigenspinors (EID, (ESD, and the unitary 
transformation ([S]) , we introduce an auxiliary matrix for 
a given value of energy at a curved surface 

Ma{0) = 

1 (im_i ) : (~^m-i)) = 

= Ui-9)MA{0). (24) 

As a result, in region I the wave function can be writ¬ 
ten as a superposition of all solutions for a curved sur¬ 
face in the form ^'/(fl) = Myi(0) exp(iTO0)C'/. Here, 
C/ is a vector of four unknown coefficients, to = 
diag(TOi, TO_i, —TOi, —TO_i) is a diagonal matrix. 

The overlap of eigenspinors of the flat and bended re¬ 
gions can be readily calculated with the aid of Fqs. m, 
(EH), which results in 

(4--)^ (^/i) ~ (-sign( 7 fc,Ao)s + a) x (25) 

X ($(()^exp(-i|CTy) ($^). 

Evidently, the overlap is zero at cr = sign{jkxEo)s. Note 
that already this result implies that some of the four 
channels between the flat and curved regions could be 
closed. 

Matching the wave functions at the boundaries of re¬ 
gions L, I, and R, for an incoming electron flow from the 


left-hand side, leads us to the following equations 
'i’LixL) = 4'7(-/)/2) ^ 

=M^(-<//2)C7, (26) 

'I'_R(a:fl) = 'I'/(-|-(/)/2) => 

Mo [ Q ] = Ma(+</)/ 2) exp(iTO0)C/ . (27) 


and r = 


We recall that the angles 9o and determine the xl,r 
coordinates: xl = Rcos{9o + (j)): xr = Rcos^q. Here, 

r(L)“ 
r(L)f 

reflection coefficients, respectively, for incoming electron 
either with a spin up | f) = [ | ) 


_ f t{L)^ 

[mi 


are transmission and 


or with a spin down 



Solutions of Fqs. (EE) (and similar equations for an in¬ 
coming electron flow from the right-hand side) yield the 
following probabilities 


|f(L)jp = |t(i?)jp = ( 28 ) 

\t{L)lf = \tiR)\f = (29) 

|r(L)|p = |r(i?)|p = l-^^^^ (30) 

|r(L)|p = |r(i?)jp = l-^^^. (31) 

Here, we have also introduced the variable 

Xy sin(TOs((i) XyR sin(TOs</) , . 

Zs — — X ,5 — ii , 

tm, 7 "is 


related to the characteristics of the curved surface (see 
Sec.HA). 

Evidently, there is no backscattering for incoming elec¬ 
trons, if Ay = 0 [see Fqs. (EE-dSE]- However, at Ay ^ 0 
backscattering with a spin inversion takes place. The 
reflection probabilities without the spin inversion are 
|r(L)||^ = |r(L)|p = 0. The same is true for the trans¬ 
mission probabilities with a spin inversion, i.e., |t(L)|p = 

|t(L):j:p = 0. Thus, backscattering with a spin inversion 
is nonzero in the ripple due to the curvature induced 
SOC produced by the Ay-term. In addition, incoming 
electrons with different spin orientations choose different 
channels (different TOs, s = ±1). 

The maximum transmission probability |t(L)|p = 1 
takes place at the condition 


m-i4>c = tto , n = 1,2,... , (33) 

[see Eos. (1281) . (|32|l ]. Evidently, this probability becomes 
dominant at the minimum transmission |t(L)|p. The 
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lowest minimum of the transmission |t(L) jp occurs at the 
condition Ea = A^;, when m+i becomes imaginary [see 
Eg. d^ ]. In other words, the propagating mode trans¬ 
forms to the evanescent mode for the channel m+i. Tak¬ 
ing into account the condition Ea = Xx in Ea. ([23|) . one 
obtains the critical angle of the curved surface (in form 
of the arc) for a maximum of spin up filter efficiency 


0c 


Trn 

TO _1 


Trny 



(34) 


where the SOC strengths Xx,y are defined by Eq. ©• 
For parameters listed in the caption of Figl2] we have 
|0c| ~ 0.9967r(n =1). For the same critical angle (pc and 
Ea = —Xx we obtain a maximum for the spin down filter 
efficiency, when m_i becomes imaginary. 

Thus, there are different channels for the spin up and 
spin down electron (hole) flows. Note that the deviation 
from the energy value Ea = could produce the equal 
transmission for spin up and spin down electrons (see 
FigEla). Therefore, it is important to choose the energy 
\Ea\ to be in the close vicinity of the energy value ~ 
Xx- For the considered parameters the filter efficiency 
is, however, small. So far this result has met with only 
limited success. 


C. Scattering model for N ripples 

To increase the efficiency we suggest connecting the 
bent parts sequentially, as shown in Fig. [TJ In particular, 
the construction with the direct-l-inverse arcs (with the 
same angle 0) transforms Eas. (I^Sl)(l?71) . to the forms 

Mo ( = Ma{-PI2)Ci , 

Mo ( Q ^ = Ma{e - 0/2) exp(-im0) x (35) 

M)^^(0/2 - 7r)M^(-|-0/2) exp(-|-im0)C'7 . 

(36) 

Since in the inverse arc the phases, accumulated from 
the point of connection with the direct arc to the point 
of connection with a straight line ( flat graphene), have 
a sign opposite that of the first one, we use exp(±im0). 

Matching the wave functions at the boundaries of re¬ 
gions L, I, II, and R, for electron coming from the left- 
hand (L) and right-hand (R) sides of the construction, 
leads us to the following nonzero probabilities 


Thus, the transmissions through one and two (di¬ 
rect-l-inverse) arcs are accompanied by the inverse 
backscattering. The considered cases imply that the 
larger the number of arcs is, the stronger the inverse 
backscattering is for one of the spin components. 

Following the recipe described in Reffl^ (interfering 
Feynman paths), and combining S-matrices for N con¬ 
nected arcs, we obtain 


= 


2 


N 


(-) 


\mi? = 


N 


N 


+c 


(-) 


N 


= \mi\\ 

) 


Here, the variable is defined as 

± Zs , s = ±l. 


(41) 

(42) 

(43) 


Evidently, at Zg = 0 the transmission probability is 1 
for any number of arcs, while Zg ^ 0 leads to the decrease 
of the transmission probability with the increase of the 
number of arcs. The suppression is, however, different for 
various transmission probabilities due to their different 
dependence on the quantum number mg. 

As shown above, conditions (l33)l . ([34|) . determine the 
dominance, in particular, of the transmission probability 
of spin up incoming electrons at Ea > 0. Indeed, a set 
Af ^ 1 of an exact replica of the consistently connected 
ripples (see FigH]) does not affect this dominance (= 1) 
for the TOg=_i channel. However, this set suppresses the 
spin down transmission probability for the rug^+i chan¬ 
nel that is proportional to a; < 1 => —>■ 0. 

We would like to point out that Eas. (|4T]) . (|^ . are valid 
for odd and even number of consistently connected rip¬ 
ples. In our model the only requirement is that the di¬ 
rect ripple has to be connected to the inverse one, the 
inverse ripple to the direct one etc. From our consider¬ 
ation it follows that, if at a certain energy, for example, 
Ea > 0 there is a high transmission probability for spin- 
up electrons from the left side of our system, one obtains 
the same magnitude for the transmission probability for 
spin-down electrons from the right side. 


III. DISCUSSION 


A. N-factor 


\t{L)\\^ = 


1 2 


l + 2(z_i)2 




|r(L)lP = 

|r(L)|p = 


l + 2(z+i)2_ 






(37) 

(38) 

(39) 

(40) 


To obtain a simple picture of the physics behind the 
enhancement of the spin filtering effect, let us consider 
the transmission at the energy Ea — Xx-, when m+i be¬ 
comes imaginary [see Eq. d^ ] and the propagating mode 
transforms to the evanescent mode for the channel m+i. 
In light of Eqs. one obtains 

R 

m+i = i — Xy = ix ,x = — ~ 0.007 (44) 

7 47 
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FIG. 3: (Color online) Dependence of transmission probabili¬ 
ties |t(L)j|^ (blue, dashed lines) and |t(I/)||^ (red, solid lines) 
on the energy Ea at fcj, = 0 for 1 (a), 20 (b), 100 (c) and 200 
(d) sequentially connected ripples (vr-arcs). The parameters 
are the same as in Fig[2l 


As a result, the variable (Ea. (l32l) ) transforms in the 
form 


z+ ~ <C 1. 


(45) 


Taking into account Eos. (14311451) . one can readily estimate 
that at iV ^ 1 


C 


(+) 

-ti 


N 


c 


(- 




-tl 




2 + {Nx(t)f => 


2 + 


(46) 

(47) 


With our choice of parameters and ^ ~ tt, this result 
yields 


|t(L)|P<l^Af> —ft!45. (48) 

The illustration of this phenomenon is displayed in Figl3] 
for the transmission probabilities through 1, 20, 100 and 
200 sequentially connected ripples (tr-arcs). Here, we 
consider the transmission as a function of the curved sur¬ 
face energy Ea of the incoming electrons (holes). A small 
difference between spin up and spin down transmission 
probabilities for one ripple (FiglS^) at Ea > 0 evolves 
to ~ 100% efficiency for the spin up transmission prob¬ 
abilities for the left-side incoming electron at N = 200 
ripples (FiglSji). The opposite picture takes place for the 
spin down transmission probabilities at Ea < 0. To re¬ 
alize such a situation one might use the Si 02 substrate 
as a gate of the curved surface, which helps control the 
concentration of charge carriers in graphene. As a result, 
one can change the charge carrier type from electron to 
holers. 


B. Spin filtering and ripple parameters R and (f) 


In light of the above analysis, without loss of gener¬ 
ality, we can consider ms4> < 1 in order to observe the 
suppression effect [see Eq. (l48l) ] . With the aid of Eq. (l23l) , 
taking into account that ^ Xy, this requirement leads 
to the following inequality 

To remain at the maximum, for example, the transmis¬ 
sion probability |<(L)jp = 1, it is necessary to fulfill con¬ 
dition (1551) . As a result, in light of Eg. (1551) and the con¬ 
dition Xx ^ Xy, taking into account Eq. dZD, one obtains 


Combining this equation with Eg. (1491) . we have 


TT — 1 , TT -b 1 

< 4 >< 


2/3 


2/3 


(51) 
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Thus, Eqs.dSOl),® determine the region of feasibility 
of the parameters R and <j), where the spin filtering ef¬ 
fects could exist at fixed system (graphene) parameters 
such as 7 , <5, and the electron energy Ea- From this ob¬ 
servation, two arguments follow in favor of our findings. 
First, even at ^ ^ (see Fa. ([34|) l one of the spin com¬ 
ponent in the incoming electron (hole) flow is suppressed 
for a large enough number of ripples at some particular 
energy region. Second, we assume that all ripples are 
identical. Practically, the graphene surface is randomly 
curved, and it is a real challenge to create identical, con¬ 
sequently connected ripples. However, it is our belief that 
modern technology will allow us to realize this situation 
soon or later. Whatever the case, the spin filtering effect 
should survive if small variations of radii and angles of 
consequently connected ripples are subject to conditions 
dsn]),®, at a fixed value Ea of the electron energy flow. 


C. Effect of a finite fcji-momentum 

In our model a single ripple is modeled as part of a 
nanotube that is infinite in the y direction. Evidently, 
realistic ripples are limited in space in both the x— and 
y— directions. In particular, graphene nanoribbons are 
considered prominent candidates to control the electronic 
properties of graphene based devices. This issue requires, 
however, a dedicated study, and is the subject of a forth¬ 
coming paper. 

In order to have some idea of what should be expected 
in graphene nanoribbons, we analyze the case with a fi¬ 
nite ky ^ 0. Nonzero ky could mimic the case of a ripple 
limited in the j/-direction. Indeed, a finite width in the 
y-direction introduces the quantization of the ky momen¬ 
tum on the curved surface. As a result, the eigenspinors 
at the curved surface would depend on the mixture Eky 
values for s = ±ljhe., altogether four momentum ky 
(see details in Refill. In this case analytical expres¬ 
sions are too cumbersome, even in a simple case of one 
conserved momentum ky on the curved surface. There¬ 
fore, we proceed with a numerical analysis that provides 
a vivid presentation of a simple case with a single value 
of the ky momentum on the curved surface. 

Let us suppose that the incoming electron flow pos¬ 
sesses a momentum k = {kx,ky) in regions L (R). Evi¬ 
dently, in this case Eq = E^^k"^ + For simplicity, 
we consider Eq > 0 , and obtain for the momentum on 
the curved surface 

fcy = <y /7 = ^^o /7 X sin(a). (52) 

The results of the calculations exhibit a degradation of 
the spin filter ability of our system. At a fixed value of 
the energy Eq = \x the transmission probability |t(L)jp 
decreases drastically at |a| > tt/S [see FigllKa)]. It seems 
that the spin filtering effects would survive at jaj < tt/S. 
Note, however, that this estimation depends on the sys¬ 
tem parameters, such as 7 and Eq. 
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FIG. 4: (Color online) Transmission probabilities |t(Z/)||^ 
(blue, dash lines) and |t(I/)|p (red, solid lines): (a) as a func¬ 
tion of the incidence angle q at Aq = Ea = Xx = 0.32eV; 
(b)as a function of the energy Ea at ky = 0.0lA“^. The 
calculations are done for 200 sequentially connected ripples 
(TT-arcs) The other parameters are the same as in Fig[2] 


At a fixed value of the momentum ky the effectivity 
of spin filtering is reduced by ^ 10% [see FigHKb)]. At 
the same time, our systems manifests a zero transmis¬ 
sion for all spin orientation for charge carriers at energies 
—0.06eE < Ea < 0.06eE due to our choice of the value 


D. The graphene purity 

We restricted our consideration to a ballistic regime. 
This approximation is well justified due to the following 
factors. The remarkable strength of the carbon honey¬ 
comb lattice makes it quite difficult to introduce any de¬ 
fects into the lattice itself. Charge impurities that could 
limit electron mobility in graphene are still an open prob¬ 
lem from both experimental and theoretical points of 
view (see, for example, discussion in ReflU). It is also well 
known that the difference in conductivity in graphene be¬ 
tween T K, Q and room temperature is no more than a few 
percent. In other words, the electron-phonon scattering 
plays a minor role. 

We recall that a typical ripple size is ~ 7 nm (see 
Refill). In our paper the ripple is modeled as the curved 








surface in the form of an arc of a circle with a radius R = 
Inm. As a result, our system length is Tri? x 200 ~ 640 
nm. Taking into account that a typical mean free path 
of electrons in single-wall nanotubes is £ ~ 1 /rm (see, 
e.g., RefH), it seems our consideration is on a reasonable 
basis. 

Thus, in our model the basic mechanism that is respon¬ 
sible for spin filtering effects is an attenuation of one of 
the transmitting modes. It transforms to the evanescent 
mode in the energy gap created by the SOC in the curved 
surface. The multiplicative action of a large enough num¬ 
ber of ripples suppresses this transmitting mode at cer¬ 
tain conditions that provide a high efficiency for the other 
one. 


IV. SUMMARY 

We have analysed the transmission and reflection of 
ballistic electron flow through a ripple in an effective mass 
approximation, when only the interaction between near¬ 
est neighbor atoms is taken into account. In our consid¬ 
eration a ripple consists of the curved surface in the form 
of an arc of a circle connected from the left-hand and 
the right-hand sides to two semi-infinite flat graphene 
sheets. Considering the curved surface as a part of the 
armchair nanotube, we have shown that the curvature 
induced spin-orbit coupling yields a backscattering [see 
Eas. (l30l31l) ] with spin inversion. This spin inversion is 
caused by the spin-orbit term that breaks spin symmetry 
(a spin projection on the symmetry axis) in the effective 
Hamiltonian of the armchair CNT. 

In the energy gap created by the curvature induced 
spin-orbit coupling there is a preference for one spin ori¬ 
entation, depending on the direction of the electron flow 
at normal incidence. The width of the energy gap de¬ 
pends in inverse proportion on the radius of the ripple. 
At this energy range the ripple acts as a semipermeable 
membrane which is more transparent for the incoming 
electrons with spin up from the left-hand side and with 
spin down from the right-hand side, and vise versa for 
the holes. In other words, there is a precursor of chiral 
transmission of spin components of the incoming electron 
(hole) flow at a fixed energy. For one ripple system this 
effect is, however, small. In order to enforce this effect, 
we extended our consideration to a curved surface of the 
sinusoidal wave type with N arcs. This step is of crucial 
importance to suppress one of the spin components and 
to support the spin inversion symmetry for the transmis¬ 
sion probability. The larger the number of consistently 
connected ripples (arcs) is, the stronger the dominance 
of a specific spin component is in comparison with the 
other in the transmission from the same direction. There 
is a cooperative effect of chiral spin transmissions pro¬ 


duced by a large number of ripples. To trace the N- 
dependence we have derived a formula for a composite 
transmission probability for well-polarized spin compo¬ 
nents: i) Eq. 6II) for spin up electrons; ii) Eg. (1421) for 
spin down electrons. Based on these results, we predict a 
strong spin filtering effect for a sufficiently large number 
of arcs in the rippled graphene system. In contrast to the 
usual waveguide that guides optical or sound waves of a 
chosen frequency in a welEdefined direction, our system 
guides spin electron (hole) waves with a well-defined po¬ 
larization in one or another direction at a certain energy. 
It seems, therefore, natural to name this system chiral 
spinguide. 

We have considered only a curved surface that owes its 
origin to an armchair nanotube. Evidently, our model 
can be extended to other types of origins. However, the 
corresponding analysis requires a separate study. We also 
neglected the effective magnetic field that arises from the 
dependence of the hopping parameter 7 on the curva¬ 
ture (see discussion in^). This effect influences the local 
density of states^!. It can cause the localization of the 
electrons on the boundary between flat graphene and the 
curved surface, similar to the boundary state for some 
types of carbon nanoparticles^^. As a result, it might 
affect the efficiency of the spinguide. Last, but not least, 
many body effects such as electron-electron interaction 
should be incorporated and analyzed as well. It is espe¬ 
cially noteworthy that electron-electron interaction, de¬ 
signed in the form of a specific potential barrier on the 
graphene sheet^^, leads to separation of spin-polarized 
states. In fact, this result is in close agreement with our 
finding, obtained for one ripple. As mentioned above, 
the curvature induced SOC simulates a penetrable bar¬ 
rier preferable for transmission of only one of two spin 
components, depending on the direction and energy of 
the incoming electron (hole) flow. It would be interest¬ 
ing to study the interplay between the SOC and electron- 
electron interaction on the electron transport in our sys¬ 
tem. Evidently, this consideration would allow us to 
study in more detail the effect of impurities on the elec¬ 
tron mobility in our system. 

In conclusion, the transparency and the mathemati¬ 
cal rigor of our results provide good grounds to believe 
that spin filtering effects found in this paper, giving rise 
to a chiral spinguide phenomenon, will be observable in 
experiment. 

Acknowledgments 

M.P. and K.N.P. are grateful for the warm hospitality 
and creative atmosphere at UIB and JINR. This work 
was supported in part by RFBR Grant 14-02-00723 and 
Slovak Grant Agency VEGA Grant No. 2/0037/13. 


^ M.I. Katsnelson, Graphene: Carbon in Two Dimensions ^ L.E.F. Foa Torres, S. Roche, and J.-C. Charlier, Introduc- 
(Cambridge University Press, Cambridge, 2012). 





9 


tion to Graphene-Based Nanomaterials (Cambridge Uni¬ 
versity Press, New York, 2014). 

® M. Katsnelson and A. Geim, Philos.Trans. R. Soc. A366, 
195 (2008) 

W. Bao, F. Miao, Z. Chen, H. Zhang, W. Jang, C. Dames, 
and C. Lau, Nature Nanotechnol. 4, 562 (2009). 

® H. Min, J. E. Hill, N. A. Sinitsyn, B. R. Sahu, L. Kleinman, 
and A. H. MacDonald, Phys. Rev. B 74, 165310 (2006). 

® T. Ando, J.Phys.Soc.Jpn. 69, 1757 (2000). 

’’’ M.V. Entin and L.I. Magarill, Phys. Rev. B 64, 085330 

( 2001 ). 

® A. De Martino, R. Egger, K. Hallberg, and C.A. Balseiro, 
Phys. Rev. Lett. 88, 206402 (2002). 

® D.Huertas-Hernando, F. Guinea, and A. Brataas, Phys. 
Rev. B 74, 155426 (2006) . 

F. Kuemmeth, S.Ilani, D.C. Ralph, and P.L. McEuen, Na¬ 
ture (London) 452, 448 (2008). 

D.V. Bulaev, B. Trauzettel, and D. Loss, Phys. Rev. B 77, 
235301 (2008). 

L. Chico, M.P. Lopez-Sancho, and M. C. Munoz, Phys. 
Rev. B 79, 235423 (2009). 

J.-S. Jeong and H.-W. Lee, Phys. Rev. B 80 075409, 
(2009). 


W. Izumida, K. Sato, and R. Saito, J. Phys. Soc. Jpn. 78, 
074707 (2009). 

M. del Valle, M. Margahska, and M. Grifoni, Phys. Rev. B 
84, 165427 (2011). 

J. Klinovaja, M. J. Schmidt, B. Braunecker, and D. Loss, 
Phys. Rev. B 84, 085452 (2011). 

K. N. Pichugin, M. Pudlak, and R.G. Nazmitdinov, Eur. 
Phys J. B 87, 124 (2014). 

M. Pudlak and R.Pincak, Eur. Phys J. B 67, 565 (2009). 
S. Datta, Electronic Transport in Mesoscopic Systems 
(Gambridge University Press, Gambridge, 2007). 

S. V. Morozov, K.S. Novoselov, and A.K. Geim, Phys. Usp. 
51 , 744 (2008). 

F. de Juan, A.Gortijo, and M.A.H.Vozmediano, Phys. Rev. 
B 76, 165409 (2007). 

T. Fujita, M.B.A. Jalil, S.G.Tan, and S. Murakami, 
J.Appl. Phys. 110 , 121301 (2011). 

L. Ying, G. Wang, L. Huang, and Y.G. Lai, Phys. Rev. B 
90 , 224301 (2014). 

A. Fasolino, J.H. Los, and M.I. Katsnelson, Nat. Mater. 6, 
858 (2007). 



